Statistical mechanics of triangulated ribbons 
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We use computer simulations and scaling arguments to investigate statistical and structural prop- 
erties of a semiflexible ribbon composed of isosceles triangles. We study two different models, one 
where the bending energy is calculated from the angles between the normal vectors of adjacent tri- 
angles, the second where the edges are viewed as semiflexible polymers so that the bending energy 
is related to the angles between the tangent vectors of next-nearest neighbor triangles. The first 
model can be solved exactly whereas the second is more involved. It was recently introduced by 
Liverpool and Golestanian Phys.Rev.Lett. 80, 405 (1998), Phys.Rev.E 62, 5488 (2000) as a model 
for double-stranded biopolymers such as DNA. Comparing observables such as the autocorrelation 
functions of the tangent vectors and the bond-director field, the probability distribution functions 
of the end-to-end distance, and the mean squared twist we confirm the existence of local twist cor- 
relation, but find no indications for other predicted features such as twist-stretch coupling, kinks, 
or oscillations in the autocorrelation function of the bond-director field. 

PACS numbers: 87.15.Aa,87.15.La,61.41.+e 
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I. INTRODUCTION 



A characteristic feature of many biopolymers is their 
high bending stiffness. Contour lengths of the order of 
fim and persistence lengths of the order of 50nm in the 
case of DNA even allow microscopy techniques to be used 
to directly observe their structure and dynamics ||[ Q. 
The model mostly used to interpret recent experimental 
data of micromcchanical manipulations of single DNA 
chains §, (§ |, |, § is that of the Kratky-Porod 
wormlike chain in which the polymer flexibility is deter- 
mined by a single length, the persistence length l p . Gen- 
eralizations account for the chain helicity and coupling 
terms between bending, stretching, and twisting allowed 
by symmetry |, [l(| Q g, || g, || || £§ |l§ [§ . 
All these continuum models of DNA neglect the double- 
stranded structure of DNA and one may ask, if this fea- 
ture could not cause qualitative different behavior. 

The bending stiffness of single- and double-stranded 
DNA, for example, differs by a factor of 25 21 . The 



simplest model which takes the double-strandedness into 
account is the railway-track model p2| where two worm- 
like chains are coupled with harmonic springs. In two 
dimensions one finds drastical consequences: the bend- 
ing fluctuations in the plane of the ribbon are strongly 
suppressed. The molecule becomes effectively stiffer on 
larger length scales. But the relevant question is: what 
are the effects in three dimensions? Liverpool et al |l], |^| 
investigated a version of the railway-track model in three 
dimensions where bending in the plane of the ribbon is 
forbidden by a constraint. Using analytical and simula- 
tion techniques they predict the existence of a low tem- 
perature regime where ribbons adopt a kink-rod struc- 
ture due to a spontaneously appearing short-range twist 



structure resulting in an oscillatory behavior of the auto- 
correlation function of the bond-director field. Further- 
more a twist-stretch coupling is predicted. 

We study the discretized version of the simulation 
model of Liverpool et al (l], ^ in the low temperature 
regime with the help of scaling arguments and MC sim- 
ulations. In order to understand and to quantify the 
effects arising from the local twist structure of the Liv- 
erpool model we compare it with an analytically more 
tractable model where the bending stiffness is defined via 
the interaction of the normal vectors so that there is no 
tendency to helical structures. Furthermore, we perform 
several MC simulation runs with an additional external 
force in order to test if the preferred buckling mechanism 
occurs via kinks. 



II. CONTINUOUS DESCRIPTION OF TWO 
COUPLED SEMIFLEXIBLE CHAINS 



A ribbon is an inextensible, unshearable rod which can 
be parameterized by the arclength s. To each point s 
one attaches a triad of unit vectors {dj(s)}. The vectors 
di(s) and d2(s) are directed along the two principle axis 
of the cross section, the vector d 3 (s) is the tangent vec- 
tor. As the triad is an orthonormal basis set they satisfy 
kinematic equations of the form 



d,:(s) = eij k Uj(s)d k (s) 



(1) 
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with eijk being the alternating tensor and Uj(s) repre- 
senting bend (wi(s) out-of-plane, and ii2(s) in-plane) and 
twist strains («3(s)) respectively. One can find a relation 
between the ordinary Frenet equations containing only 
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two parameters, the curvature k(s) and the torsion r(s) 



dt(s) 
ds 

db(s) 
ds 

dn(s) 

ds 



K(s)n(s) 
-r(s)n(s) 
r(s)b(s) - «(s)t(s) 



(2) 
(3) 
(4) 
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and Eqs. (Q) by fixing d 3 (s) = t(s) so that di(s) and 
d 2 (s) are given by a rotation around t(s) with angle 
^(s). In this context ^(s) can be seen as the twist 
angle |l7], B3|. A straightforward calculation gives for 
the generalized torsions: u\{s) — — ^d 3 (s) • d 2 (s) = 
k(s)cos*(s), u 2 (s) = £d 3 (s) ■ di(s) = k(s) sin *(s), 

and ua( a ) = £di(«) . d 2 (s) = r(s) + The total 

Twist Tto of a ribbon is thus given by the integration of 
the local twist us(s) along the contour normalized by the 
factor 2ir 



Tw = — [ u 3 (s)ds 
2tt Jo 



(5) 



with L being the contour length. Together with the pa- 
rameter set Ui(s), which determines whether the stress- 
free reference configuration includes spontaneous curva- 
ture and twist, the elastic part of the Hamiltonian is usu- 
ally defined by quadratic terms in Ui(s) — iii(s) |l2] , 

HO, 0, H H H, H, 

It is an interesting question to which extent this 
generic description applies to more microscopic models of 
DNA (25) . The simplest case is that of a "railway track" 
or ladder model consisting of two (or more) semiflexible 
chains 



Hi 



ds 



ds 2 



d 2 T 2 (s y 

ds 2 



(6) 



plus a coupling between opposite points on different 
chains p^j. Liverpool et al pi || considered the limit 
where the distance a between the coupling points (i.e. 
the width of the ribbon) is imposed as a rigid con- 
straint which prevents bending in the plane of the rib- 
bon: ^jjp • b(s) = where t(s) = ^j^- is the tangent 
vector to the mid-curve r(s) = ri(s) — | = r 2 (s) + | 
and b(s) is the bond-director pointing from one strand 
to the other. Note, that the constraint is equivalent to 
^(s) = 0. Rewriting Eq. (j^) in terms of ribbon variables 
they found 



Hi 



ds < 2 



d 2 r(s) 
~d^~ 



a 

y 



d 2 b( s y 

ds 2 



(J) 



FIG. 1: Illustration of the used variables. The length of each 
triangle |ti| corresponds to the bond length b and the height 
bi| = i&tan(a) defines the strand separation length. {6i} 
term the folding angles. 



which can also be expressed as 



/dt 

Us 



d 2 h 



du\ 
ds 

dn\ ' 
ds) 
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(n 2 -r 2 f. 



(8) 



(9) 



Note, that henceforth we use b(s) as the bond-director 
and n(s) as the normal vector to the ribbon plane. 



III. GEOMETRY OF TRIANGULATED 
RIBBONS 



Following Liverpool et al |l], ^ we consider ribbons 
discretized by triangulation. In order to extract some 
fundamental properties of double-stranded semiflexible 
polymers we consider a ribbon-like system composed of 
isosceles triangles as shown in Fig. [lj The orientation 
of each triangle is given by N — 1 rotations around the 
edges of the triangles with folding angles {#;}. N is the 
number of triangles characterized by a set of trihedrons 
{ti, bi, rii} where ti is the tangent vector of the ith trian- 
gle, bi is the bond-director, and rii is the normal vector. 
Going from one set of trihedrons {ti, bi, rii} to the neigh- 
bor set {ti+i,bi+i, ni+i} implies a rotation TZi around 
the edge between the respective triangles with angle 0i 
and a reflection of bi and rii, he. 




with 



T 




ti • ti + i ti • bi + i ti • rii+i 
bi • t i+ i b ; • b i+ i bi • n i+ i 
rii -ti+i rii -bi+i rii -rii+i 



(10) 

(11) 
(12) 
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The matrix product TIZi can be viewed as a transfer 
matrix. The evaluation of the scalar products of IZi gives 

TZi,U = cos(0j) + cos(a) 2 (1 — cos(0j)) 

TZi t %2 = — cos(a) sin(a) (1 — cos(0i)) 

72-i,i3 = - sin(a) sin(0i) 

72-i,2i = cos(a) sin(a) (1 — cos(0j)) 

TZ il2 2 = cos(^) + sin(a) 2 (1 - cos(0;)) (13) 

72 4 ,23 = - cos(a) sin(6»j) 

72-1,31 = sin(a) sin(6>;) 

72^32 = cos(a) sm(6i) 

72j,33 = cos(^). 

In order to quantify properties such as bending and 
twisting within the given discretization we study the re- 
lation between the folding angles Oi and these quantities 
which is illustrated in Fig. |^. One recognizes that the 
chain is not bent in case of B\ — flj+i = 50 i = and that 
one gains purely twisted structures if Oi = const. On the 
other hand the chain is untwisted but bent if 50 i = 28 \. 
In case of Oi ^ ±0j+i and 6*^=^0 the chain is bent and 
twisted simultaneously resulting in solenoidal/torsional 
structures as is illustrated in Fig. ||(f). A kink is charac- 
terized by unlike twists meeting at an edge as it is shown 
in Fig. |(d). 

Due to the triangulation of the ribbon one has to 
consider three triangles to calculate the discrctizcd ex- 
pressions for the out-of-plane bending strain u±(s) — 
-£t(s).n(s) « _ t(^)-t(,) , n(g) = __i_ t(s+As) . n(s) 

and the twist strain u 3 (s) — 4-b(s)-n(s) » b (s+As)-b(s) _ 
n(s) = ^b(s + As) • n(s) which we call Ki and Tj respec- 
tively. The local curvature Ki and the local twist rate t% 
between triangle i and i + 2 are therefore given by 

Ki = -- }^ nj ■ t j+1 « —j^-Mi (14) 
n = I^ nj . bj+1 «H2^( ( 9. + ( 9<+1 ) ) (is) 

where the accuracy of the right-hand side expressions 
only depends on the refinement of the discretization, i.e. 
on the values of b and a. Hence a spontaneous bending 
can be introduced via an additional term to the Hamil- 

tonian with H CU rv = k curv ^ (j2]=\ n j ' t j+i - S0 sp ^ 
and a spontaneous twist can be introduced by an ad- 
ditional term Ti Tw = k Tw £\ fejii n i ' b i+i _ Qs P ,i\ ■ 
Note, that the total twist Tw is given by Tw — 
1/(271") E^. 

IV. MODEL DESCRIPTION 

The bending stiffness within the given discretization 
can be taken into account by various interactions. One 




(a) (b) (c) (f) 



FIG. 2: Illustration of bending, twisting, and kinking, (a) A 
flat ribbon as ground state conformation, (b) A twisted struc- 
ture (c) The same twisted structure obtained with a smoother 
discretization, (d) Unlike twists meeting at the center re- 
sulting in a kink with 8i positive for i < N/2, negative for 
i > N/2, and \6i\ = |0 i+1 |, i.e. 56i = 0, Vi ^ N/2 and 
58n/2 = 20JV/2- (e) A bent structure, (f) A mixture of bent 
and twist resembling a solenoidal structure. 

possible definition of a bending stiffness, which makes 
the problem analytically tractable, is a nearest neighbor 
interaction (plaquette stiffness) between the normal vec- 
tors {rii} in analogy to the triangulation of vesicles |2(| 
which results in the following Hamiltonian 

nj N-l 

^ = fc£(l + ni .n i+1 ). (16) 

^ i—l 

In contrast Liverpool et al jl[ |) were interested in the sta- 
tistical mechanics of coupled wormlike chains and there- 
fore chose a next-nearest neighbor interaction (edge stiff- 
ness) between the tangent vectors {ti} with rigidity k so 
that the Hamiltonian is given by 

= *E(l-tr*H. a ). (17) 

i—l 
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Both definitions lead to a flat ribbon as the ground state 
conformation for zero temperatures T = 0. 

The above defined interactions lead to very distinct 
conformational features of the ribbon which can be un- 
derstood by building up the ribbon just by adding suc- 
cessively the triangles in the absence of thermal fluctua- 
tions. Assuming that 9\ ^ all subsequent angles 9i with 
i > 1 vanish in the case of the nearest neighbor interac- 
tion (H nn ). In contrast the tangent-tangent interaction 
(Htt) leads to the formation of a helix with 9i = 9 i+ i 
as a result of the enforced alignment of the tangent vec- 
tors. This suggests a correlation of the folding angles 
{9i} which entails at least locally helical structures. 

Assuming that the chains are rather stiff (continuum 
limit), i.e. small folding angles one can expand the 
Hamiltonians with regard to 6i. Since H nn is diagonal in 
9 il it is sufficient to consider terms up to second order. 
Htt contains coupling terms between 9i and 9i + i which 
makes it necessary to keep terms up to fourth order in 
the analysis: 



N-l 



k B T 
Htt 



i=l 

i=i ^ 
+ sin(a) 2 cos(a)X 2 ^ +i } 



(18) 



(19) 



with 



MC SIMULATION 



Both models have local interactions and can be studied 
conveniently using a dynamic MC scheme. Trial moves 
consist of small random changes of the folding angles by 
a small amplitude 1/ Vk : where k is the bending stiffness, 
and are accepted or rejected according to the Metropo- 
lis scheme |27j . In the simulations we always use the full 
Hamiltonians Eq. @ and |[7|). MC moves changing the 
folding angles correspond to the well-known Pivot algo- 
rithm 1 28 . The conformations are subsequently recalcu- 
lated from Eqs. (10)-(O|) and analyzed. Each simulation 



run comprises 100000 MC-moves where one MC move 
corresponds to N — 1 trials with N being the number of 
triangles. The longest correlation time we observed was 
of the order of 50 MC moves for the total twist of the 
chain. In order to check if equilibrium is reached we com- 
pared simulation runs with a flat initial conformation, i.e 
9i = 0, with simulation runs with crumpled conforma- 
tions corresponding to equally distributed angles 9i out 
of the interval [—1/Vk; 1/y/k]. Both runs yield the same 
results for the calculated observables. 



VI. PLAQUETTE STIFFNESS 

Since the Hamiltonian H nn of Eq. ( |l8| ) is quadratic 
and diagonal in 9i the solution in angle space is triv- 
ial. As a consequence of the independence of succes- 
sive folding angles it yields (9i9j) — jSij and (A) = 

(Ui=i( Tn k)) = (TTZk)^ 1 where the matrix product is 
carried out in the eigenvector basis of (TTZk) (the eigen- 
vectors depend only on the geometry of the triangles). 
The diagonal elements of (A) are the correlation func- 



tions of (ti • tj), (bi • bj), (rij • rij). Thus one calculates 



(TIZk), diagonalizes it, raises it to the power of j — i, 
transforms it back, and performs the continuum chain 
limit with s = (j — i)b, l p = bkj sin(a) 2 , a = i&tan(a), 
U~ *) ~~ * 00 j b — > 0, i.e. a — > 0, where l p is the persistence 
length, a is the strand separation, b is the Kuhn segment 
length, < s < L is the arclength, and L is the contour 
length. Note that within this model a is a fixed parame- 
ter that determines bending characteristics of the ribbon. 
An exact expression for the autocorrelation functions is 
obtained: 



<t(o).t(»)> 

(b(0)-b(,s)) 
(n(O)-nOO) 



exp 
exp 
exp 



l p tan(a)' 
s 

l p sin(a) 2 



(20) 
(21) 
(22) 



For a = 7r/2 one recovers the usual wormlike chain result 
for two dimensions. All crosscorrelation functions (the 
off-diagonal elements of (A)) vanish. Eqs. ( f2ll) , ( |22| ) rep- 
resent the persistence length l p ^ n — l p tan(ap for bend- 
ing within the plane of the ribbon and the persistence 
length l p , ou t = l p sin(a) 2 for bending out of the plane 
of the ribbon respectively 0] . This model was recently 
treated as a twisted zig-zag fiber within the framework 
of a two-angle model for studying structural properties 



of chromatin |29 



From the tangent-tangent correlation function one can 
calculate the mean squared end-to-end distance: 

R% = <(R(L) - R(0)) 2 ) = / dsi / d* 2 <t(si) • t(« 2 )) 

Jo Jo 



■211 „ ~ -111 j I - exp j -y 



(23) 



Eqs. ( pp| ) and (^3|) are identical to results for single worm- 
like chains [3C| ]. Eq. ( p3| ) interpolates between the limit- 
ing behaviors of random coils (2Ll p ) for L> l p and rigid 
rods (L 2 ) for L <C l p . 



VII. EDGE STIFFNESS 

In the following we present a simple scaling argument 
which allows us to rationalize the behavior of the Liver- 
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FIG. 3: (a), (b) Scaling plots for (0?) (downward triangles), 
(59f) (squares), and g with a — 7r/3 (upward triangles) and 
k — 50 respectively, (c) Numerical evidence for the derived 
expression of (0i0j). The data refer to k = 50 (squares), 
k — 100 (circles), k = 200 (upward triangles), k — 500 (down- 
ward triangles) and a = n/4. We determine the correct pref- 
actor gfu = 0.56 ± 0.05 of g from the numerical data of the 
folding angle correlation function (8i0j) which is our solely 
free parameter and use it for all following comparisons be- 
tween scaling analysis and numerical results. 



pool model. Consider first the 89 part of Eq. (19). In the 



absence of other terms the folding angles would perform 
a simple random walk with step length (66?) = ksi ^ a ^ ■ 
The leading term limiting the fluctuations of the fold- 
ing angles around zero is of order O(0f). The behavior 
of the coupled system can be inferred from scaling ar- 



guments similar to those used for polymer adsorption. 
Consider a vanishing folding angle and follow the chain 
in either direction. Up to a characteristic number of steps 
g the folding angles will show simple diffusion. As a con- 
sequence the mean-squared folding angle averaged over 
this short segment is (6?) = g(89f) corresponding to a po- 
tential energy |- ~ g{0f) ~ 3<?(0f) 2 ~ 3 5 3 (^|). The 
free diffusion of the folding angles has to stop when this 
potential energy is of order fc^T, suggesting 



1 



k sin(a) 2 
k tan(a) 2 



exp 



2.9 



(24) 

(25) 
(26) 
(27) 



Fig. H shows that these arguments are fully supported 
by the results of our MC simulations with g = (0.56 ± 







.05) ( 



k tan(a)' 
3 



Using again the transfer matrix ansatz one obtains in 
the low temperature limit by considering only terms of 
the order 0(9f) the following expression for ti - tj , bi • bj 
and rij • nj: 



ti • tj 



bi-bj = 1 




(28) 



/ n2 /J - 1 J -1 3—1 



\k—i 



J/2 



k=i k'=k+l 



3-1 3-1 



(29) 



1 



k=i/2 



k—i k'—k-\-l 



(30) 



Note that i,j are either odd or even depending on which 
strand is under consideration. Without loss of general- 
ity we choose i,j to be even. First of all we use that 
(t(0) • t(s)) has to interpolate between 1 for s = and 
for s — > oo and that the right hand side of Eq. d2^) 
is the Taylor expansion up to first order of the exponen- 
tial function exp ^ hin ^ a ) ^J^ 2 ^ <5#2fc) ^ ■ Substituting 

then s = 2\j — i\b and l p = Abk, performing the contin- 
uum chain limit with 6^0 and a — > 7r/2 respectively, 
i.e. keeping the strand separation a constant, yields the 
following expression for the autocorrelation function of 
the tangent vectors: 



(t(0)-t( S ))=exp ( -- 



(31) 
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FIG. 4: Comparison of MC-data and analytical results (solid 
line) for the autocorrelation function of the tangent vectors 
with k — 50 (squares), k — 100 (circles), k — 200 (upward 
triangles), k — 500 (downward triangles) and a = 7r/3, and 
k — 50 and a = tt/4 (diamonds), a = tt/6 (pentagons). 



Thus the mean squared end-to-end distance R 2 E becomes 
identical to Eq. (p3|). Eq. (|3l|) is confirmed by our MC 
simulation data shown in Fig. |j. 

To get an idea of the structural properties character- 
ized by the autocorrelation function of the bond-directors 
(bi • bj) we calculate the mean squared twist (Tw(i,j) 2 } 
of the ribbon. Following the definition of the local twist 
rate r, of Eq. ( |i~5| ) the total twist between two triangles 
of index i and j is just the sum of the local twist angles 
determined by the projections of the normal vector of 
the ith triangle onto the bond-director of the (i + l)th 
triangle, that is 



1 

Tw(i,j) = - Vrii -b i+ i 

Z7T » — ' 



cos(a) 
2tt 



3-1 



k ■ 



(32) 



Comparing Eq. (29) and ( p2|) we find for small twist an- 
gles 



(b i -b i ) = l-2n 2 (Tw(i,j) 2 ). 



(33) 



Hence the autocorrelation function of the bond-directors 
can be seen as a measure for the local twist structure of 
the ribbon. 

In contrast to the plaquette stiffness model the an- 
gles 9i are correlated (see Eq. (P7|)). Therefore the dou- 
ble summation over {6iQj) in Eq. ( p9| ) yields an anal- 
ogous result as it is obtained in the calculation of the 
mean squared end-to-end distance of the wormlike chain 
model. Using the derived scaling expressions of Eqs. ( |25| ) 
and (|26|), the same substitutions as in Eq. (|3l|), and per- 
forming the continuum chain limit one obtains the fol- 
lowing relationship for the autocorrelation function of the 
bond-directors: 



(b(0) ■ b(s)) = exp (-2vr 2 (Tu;(0, sf)) 



(34) 




FIG. 5: (a) Autocorrelation function of the bond-directors 
with k = 50 (squares), k — 100 (circles), k = 200 (upward tri- 
angles), k — 500 (downward triangles) and a = tt/3. The data 
show the predicted functional (solid line) form for (b(0) -b(s}) 
of Eq. rt34|) . In order to check the scaling argument of Eq. (134) 
we determined the correct prefactor gfu = 0.56 ± 0.05 of g 
with the help of the numerical data of (Oidj) (see Fig. ^ and 
put it in Eq. (|34|). The agreement is excellent. The dashed 
line which oscillates is the predicted functional form of Liver- 
pool et al §, §. 

(b) Comparison of our simulation data with the analytical 
result of Liverpool et al (dashed line). The predicted oscil- 
lation and resultant pitch is not recovered. But we find the 
same scaling behavior of the helical persistence length with 

i 2 

l b = gb ~ lp a 3 . It is also striking that the predicted func- 
tional form of Liverpool et al is in very good agreement with 
our numerical data within one helical persistence length Z&. 



with 

(Tw(0,s) 2 ) 



37r 2 \gb 



2 1 



exp 



s 

'2g~b 



(35) 



and gb — g l J^lp ai /3 1 / 3 , where gfu = 0.56 ± 0.05 is the 
fitted prefactor for the scaling function g. a represents 
the strand separation of the ribbon which is given by 
a = |bi| = ifotan(a). Hence we observe two length scales 
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FIG. 6: Autocorrelation function of the normal vectors with 
k — 50 (squares), k — 100 (circles), k — 200 (upward tri- 
angles), k — 500 (downward triangles) and a = it/3. We 
divided out of the normal vector correlation function (solid 
line) the tangent correlation function (t(O)-t(s)} (see Eq. (^6|)) 
so that one should regain the same exponential decay as for 
(b(0) ■ b(s)} which is in agreement with the numerical data. 



influencing the local twist structure of the ribbon: on the 
one hand the single strand persistence length l p and on 
the other hand the strand separation a. The predicted 
scaling behavior of (b(0) • b(s)) can be observed in the 
simulation data as it is shown in Fig (||). Note that (b(0)- 
b(s)} as well as all other calculated observables within 
this model is independent of the geometry of the triangles 
in contrast to the previous model where a was a fixed 
parameter. 

Eq. (|30| ) can be evaluated in the same manner resulting 

in 



(n(0)-n( S )) = (t(0)-t( S ))(b(0)-b( S )) 



exp 



2ir 2 (Tw(0,s) 



(36) 



Eq.(p4[) shows that the autocorrelation function of the 
normal vectors is the product of (t(0) • t(s)) and (b(0) • 
b(s)). For very stiff chains, the tangent correlation func- 
tion gives just small corrections to the normal vector cor- 
relation function. Therefore one can interpret Eq. (|34|) as 
the rigid rod limit of Eq. (|3^) . 

Other important structural properties of the ribbon 
can be extracted out of the crosscorrelation functions. 
(n(0) -t(s)) and (n(0) -b(s)) describe the mean curvature 
and mean twist respectively and vanish in both models 
for symmetry reasons. For (b(0) • t(s)) we empirically 
observe the following relationship: 

(b(Q)-t(s)) = (2^) 2 a ^(T^0, S )V (2 " )2<7M °" )2> - 
as 

(37) 

Eq. ([37j) can be understood qualitatively in the following 
way. Due to the anisotropic rigidity of the ribbon the 



FIG. 7: Crosscorrelation function of the bond-directors and 
the tangent vectors with k — 50 (squares), k — 100 (circles), 
k — 200 (upward triangles), k — 500 (downward triangles), 
k — 1000 (diamonds) and a = n/3. The data validate the pre- 
dicted functional form (solid line) for (b(0) -t(s)} of Eq. 



scalar product b(0) ■ t(s) is only non-zero if the chain 
is bent and twisted simultaneously. In case the ribbon 
is either solely bent or solely twisted the bond-directors 
are always perpendicular to the tangent vectors and the 
scalar product b(0) • t(s) vanishes for all s. The rate of 
mean twist of one helical persistence length lb = gb which 
defines the size of the locally existing helical structures 
can be calculated with Eq. at 



_ yielding y/(Tw(0, bg) 2 ) w 
ijij. This corresponds to a typical twist angle of W = -| 
using Tw — 27r>F . Within 4 the twist rate is determined 
by the derivative of the mean squared twist j^(Tw(0, s) 2 ) 
which gives rise to the increasing correlation function 
(b(0) • t(s)) up to the maximum value at If, = gb. For 
larger internal distances of the chain the rate of mean 
twist is a random sequence of ± so that the crosscorre- 
lation function has to vanish and therefore decreases ex- 
ponentially with exp (— (27r) 2 (Tw(0, s) 2 )) . Fig. (0) com- 
pares Eq. ([37]) with our numerical data. The agreement 
is excellent. 



VIII. 



BEHAVIOR UNDER COMPRESSION: 
EULER BUCKLING VS. KINKS 



As discussed in section III the edge stiffness model in- 



cludes local twist correlations at least on small length 
scales as a consequence of the correlation of the folding 
angles {0i}. In order to understand and to quantify the 
effects arising from the local twist we measured the prob- 
ability distribution functions of the folding angles, of the 
twist, and of the end-to-end-distance for different rigidi- 
ties and compared the latter with the usual wormlike 
chain model to see which differences occur. 

If there is a preference for kinking one can enforce 
this property by applying an additional constant force 
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FIG. 8: Comparison of the probability distribution functions 
of the total twist of the ribbon for both models with / = 
/ = {0, 0.01, 0.02, 0.09} and l p = L = 400 with the scaling 
analysis for / = 0. One recovers the same Gaussian shape for 
all values of /. 



Fbuck = /Rb/-Rb which compresses the ribbon. In ad- 
dition the change in the end-to-end distance Re caused 
by the buckling force should affect the twist distribution 
function P(Tw) if Re and Tw are coupled. 

For small forces we calculate the change of twist under 
the influence of the external force F^uck — f within the 
framework of linear response theory: 

(ATw(0,L) 2 ) = (Tw(0,L) 2 ) - (Tw(0, L) 2 ) f=0 

= -13 f ((R E Tw(0, L) 2 ) f=0 - (R E ) f=o(Tw(0, L) 2 ) f=0 ) 

(38) 

with P = X/ksT. This predicts a change of 
the mean squared twist of the chain if a twist- 
stretch coupling determined by (ReTw(0, L) 2 ) f =0 ex- 
ists. Note that (ReTw(0, L))/=o vanishes due to 
symmetry reasons. The evaluation of our numeri- 
cal data yields that (ReTw(0,L) 2 ) is uncorrelated, 
too. To quantify if higher order terms in / con- 
tribute to a change of (Tw(0, L) 2 ) we carried out 
several simulation runs with varying force strengths 
/ = {0,0.01, 0.02,. ..,0.09} corresponding to ^§ = 
{1, 0.95, 0.87, 0.71, 0.51, 0.36, 0.26, 0.21, 0.17, 0.15}^ 

Fig. | shows the same Gaussian shape for all measured 
probability distribution functions of the total twist of the 
ribbon P(Tw(0, L), /). This implicates that there is no 
twist-stretch coupling inherent in the system. The same 
is valid for the distribution function of the folding angles. 

Moreover we measured the probability distribution 
function P(Re, /) of the end-to-end distance Re for all 
applied forces /. Using the multiple histogram method 
developed by Ferrenberg and Swendsen |[l] one can then 
recombine all measured histograms with a reweighting 
procedure to a single probability distribution function 
P(Re) with overall very good statistics. Fig. |^ shows 



FIG. 9: Probability distribution functions of the end-to-end 
distance of the edge stiffness model for different discretiza- 
tions (N = 800 upward triangles, TV = 600 downward trian- 
gles) with a = — and l p = L — 400 calculated with the help 
of the multiple histogram method and the usual worm- 
like chain model (solid line) . The PDF of the wormlike chain 
model is calculated with the derived analytical expression of 
Wilhelm and Frey 53. 



P(Re) for Ti-tt and the wormlike chain model. Quite 
contrary to a shift to noticeably shorter end-to-end dis- 
tances Re as one would expect for the above described 
phenomena of kinks one just recovers the usual worm- 
like chain behavior. This indicates that the ribbon just 
bends under the external force in contradiction to a kink- 
rod structure. Another quantity which is sensitive to 
the presence of kinks is a three-point correlation function 
of the end-to-end distance Re and the twist to the left 
Tw(0, -§), and to the right Tw(^,L) of the center of the 
chain. Due to the buckling force the center of the chain is 
labeled which means that a kink is detected if the end-to- 
end distances with Tw(0, ^)Tw(^,L) < (unlike twists 
meeting at the center) are smaller than the end-to-end 
distances with Tw(0, j)Tw(^,L) > (like twists meet- 
ing at the center). Fig. [lO] shows the mean end-to-end 
distance depending on the value of Tw(0, ^)Tto(|, L) for 
l p = 200, L = 400, and / = 0, / = 0.03, / = 0.06. We 
do not find an asymmetry between like and unlike twists 
meeting at the center as it would support the prediction 
of kinks made by Liverpool et al (3, 



IX. SUMMARY 

We have reinvestigated the mechanical properties of 
the model introduced by Liverpool et al 0, | of a double- 
stranded semiflexible polymer and rationalized the re- 
sults of our MC simulations with the help of a simple 
scaling argument. We recover the predicted simple ex- 
ponential decay of the tangent-tangent correlation func- 
tion with the single strand persistence length l p and that 
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FIG. 10: End-to-end distance Re as a function of the prod- 
uct of the twist left and right of the center of the chain 
Tw(0,%)Tw(%,L), which is a measure for unlike (nega- 
tive sign) and like (positive sign) twists meeting at the 
center, and as a function of the applied buckling force 
with / = (squares), 0.03 (circles), 0.06 (triangles) and l p = 
0.5L = 200. Re refers to the average of one interval of 
Tw(0, ~)2\u(§ ,L) and (R E ) refers to the mean value of all 
sampled end-to-end distances. One does not find an asym- 
metry between end-to-end distances for like and unlike twists 
meeting at the center. The larger fluctuations for larger val- 
ues of Tw(0, L)Tw{^,L) are the result of a poorer sampling 
rate. 

t(0) • t(s) is independent of the separation a of the two 
strands, which is in addition to l p the other relevant 
length scale in the problem. Also in agreement with 
Ref. jl], |^] we find that the helical persistence length 4 
and the helical pitch P scale with l^ 3 a 2 / 3 . Qualitatively, 
one would expect to see oscillations in the bond-director 
correlation function, if P < If,. This can be understood 
by calculating the rate of mean twist within lb = gb, i.e. 
\J (Tw(0,gb) 2 ). If the mean twist rate exceeds n an oscil- 



latory behavior has to be observed. But our calculation 
gives a twist rate within lb = gb of approximately ±1/16. 
For larger distances of the chain the rate of mean twist is 
just given by a random sequence of ±1/16 and thus can- 
not account for an oscillatory behavior of (b(0) • b(s)). 
Liverpool et al predict P = l bj while our analysis indi- 
cates P = 16h as it is demonstrated in Fig. || (b). The 
authors claimed support from their own simulations, but 
failed to provide a quantitative comparison between their 
numerical and analytical results. In fact the presented 
oscillations seem to be ordinary fluctuations within the 
statistical errors. But as can be seen in Fig. | (b) the 
predicted functional form for the bond-director autocor- 
relation function is in very good agreement with our nu- 
merical data as well as with our scaling results within 
one helical persistence length 4 = gb. 

Moreover our simulation results with applied con- 
stant buckling forces do not provide any evidence 
of the predicted tendency of kinking or the claimed 
twist-stretch coupling. Thus contrary to the claim 
made in Ref. [[[J |) the local twist structure does not 
suffice to explain experimental observations such as 
the twist-stretch coupling || Q and the kink-rod 
structures (34| of helical double-stranded molecules. 
These features require the inclusion of a spontaneous 
twist incorporated by an additional term in the Hamil- 

tonian, e.g. H Tw = k Tw J2i fejii n j 1 b i+i _ ^s P ,i) i 
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